function [PV_Mean_Dist, Shuffle_PV_Mean_Dist] = PV_Mean_Dist_One_Animal_One_Dose(Data_Summary, drug, animalID)

for n=1:size(Data_Summary,2)
    if contains(Data_Summary(n).drug, drug) && contains(Data_Summary(n).animalName, animalID)
        PV_Dist(n,:)=Data_Summary(n).PV_Proportion;
    end
end
for n=1:size(PV_Dist,1)
    if max(PV_Dist(n,:))==0
        PV_Dist(n,:)=nan;
    end
end
for n = 1:size(PV_Dist,2)
    PV_Mean_Dist(1,n) = mean(PV_Dist(:,n), 'omitnan');
end


for n=1:size(Data_Summary,2)
    if contains(Data_Summary(n).drug, drug) && contains(Data_Summary(n).animalName, animalID)
        Shuffle_PV_Dist(n,:)=Data_Summary(n).Shuffled_PV_Proportion;
    end
end
for n=1:size(Shuffle_PV_Dist,1)
    if max(Shuffle_PV_Dist(n,:))==0
       Shuffle_PV_Dist(n,:)=nan;
    end
end
for n = 1:size(PV_Dist,2)
    Shuffle_PV_Mean_Dist(1,n) = mean(Shuffle_PV_Dist(:,n), 'omitnan');
end

end